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ABSTRACT 

We report on a 350-ks NuSTAR observation of the magnetar IE 1841—045 taken in 2013 September. 
During the observation, NuSTAR detected six bursts of short duration, with Xgo 5 1s. An elevated 
level of emission tail is detected after the brightest burst, persisting for ~lks. The emission showed 
a power-law decay with a temporal index of 0.5 before returning to the persistent emission level. The 
long observation also provided detailed phase-resolved spectra of the persistent X-ray emission of the 
source. By comparing the persistent spectrum with that previously reported, we find that the source 
hard-band emission has been stable over approximately 10 years. The persistent hard X-ray emission 
is well fitted by a coronal outflow model, where pairs in the magnetosphere upscatter thermal 
X-rays. Our fit of phase-resolved spectra allowed us to estimate the angle between the rotational and 
magnetic dipole axes of the magnetar, a mag = 0.25, the twisted magnetic flux, 2.5 x 10 26 G cm 2 , and 
the power released in the twisted magnetosphere, Lj = 6 x 10 36 erg s _1 . Assuming this model for 
the hard X-ray spectrum, the soft X-ray component is well fit by a two-blackbody model, with the 
hotter blackbody consistent with the footprint of the twisted magnetic field lines on the star. We also 
report on the 3-year Swift monitoring observations obtained since 2011 July. The soft X-ray spectrum 
remained stable during this period, and the timing behavior was noisy, with large timing residuals. 

Subject headings: pulsars: individual (IE 1841—045) - stars: magnetars - stars: neutron - X-ray: 
bursts 


1. INTRODUCTION 

Magnetars are neutron stars to have emission that 
is po wered b y the decay of th eir intense magnetic 
f ields (iDuncan fe Thompsonl[l992i ; iThomnson fc Duncanl 
1199611 . The magnetic field strengths inferred from the 
spin- down parameters are typically greater than 10 14 G 
(e.g., iVasisht fe GottheIlll997t iKouveliotou et all 1 1 9981) . 
though there are several sources with lower inferred field 
strengths (e.g., SGR 0418+5729, Swift J1822 3-1606: 
IRea et all 120101 : iScholz Kaspi fc Cummingl 12014 (1. There 
are 28 mag netars discovered to dat e, including six can¬ 
didates (see lOlausen &; KasT)ill20 l lTl fl 
Short X-ray bursts are often detected from magnetars. 
The bursts have a variety of morphologies, spectra, and 


1 See the online magnetar catalog for an up- 
to-date compilation of known magnetar properties, 
ht t p: / / www. physics. mcgill. ca / ^pulsar/magnetar/main. ht ml 


energies, and are though t to be produced by cr ustal or 
magnetospheric activity (IThomnson et al] 12002 (1. Inter¬ 
estingly, some bursts are followed by a long emission tail 
while others are not. It has been suggested that the en¬ 
ergy in the burst and the integrated energy in the burst 
tail are correlated (for SGR 1900+14 a nd SGR 1806—20, 
iLenters et all 120031 : iGbgfis et all 1201 ll) . which may im¬ 
ply tha t their relative stre ngths show a narrow distri¬ 
bution. IWoods et al. (]2005l ) suggested a bimodal distri¬ 
bution for the relative strengths across magnetars, im¬ 
plying two distinct physical mechanisms for the bursts. 
In addition to short X-ray bursts, giant flares and dra¬ 
matic increases in the persistent emission over days to 
months have been also observed in some sources (see 
Woods fc Thom psonll2006l : iKasni 1120011 : IRea fc Esnositol 

20111 : lMereghcttill2013l . for reviews). 

The persistent emission of magnetars in the X- 
ray band below ^10 keV is dominated by thermal 
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Table 1 

Summary of observations used in this work 


Observatory 

Obs. ID 

Obs. Date 
(MJD) 

Exposure 

(ks) 

Mode a 

Chandra 

730 

51754.3 

10.5 

cc 

XMM-Newton 

0013340101 

52552.2 

3.9 

FW/LW 

XMM-Newton 

0013340102 

52554.2 

4.4 

FW/LW 

Chandra 

6732 

53946.4 

24.9 

TE 

Swift h 

00080220003 

56241.3 

17.9 

PC 

Swift h 

00031863050 

56551.3 

4.3 

WT 

Swift h 

00080220004 

56556.5 

1.9 

PC 

NuSTAR 

30001025002 

56240.9 

48.6 


NuSTAR 

30001025004 

56540.3 

35.9 


NuSTAR 

30001025006 

56542.4 

77.9 


NuSTAR 

30001025008 

56547.1 

85.7 


NuSTAR 

30001025010 

56549.6 

53.3 


NuSTAR 

30001025012 

56556.5 

100.7 



a CC: Continuous Clocking, FW: Full Window, LW: Large Window, 
PC: Photon Counting, TE: Timed Exposure, WT: Windowed Tim¬ 
ing. MOSl,2/PN for XMM-Newton. 

*'Swift observations used for spectral analysis in Section 13.2.31 Re¬ 
sults of analysis on a larger Swift monitoring dataset are presented 
in Section |4] 

2. OBSERVATIONS 


emission and is often modeled with two blackbod- 
ies from two hot regions on the surface, or with a 
surface blackbody plus a po wer law resulting from 
magnetospheric reprocessing dThompson ct all 120021 : 
iZane et all 120091) . Some magnetars also show signifi- 
cant emission in the hard X-ray band above ~10 keV 
dKuiper et al.l 120061) which i s believed to be produced 


in 

the magnetosphere ([Thompson & Beloborodov 

2005 

Hevl & Hernauistl 20051: Baring & Harding 

2007 

: Beloborodov & Thompson! 120071). Recently, 

Beloborodov! (120131) proposed a coronal outflow model for 


the hard X-ray emission. The model makes specific pre¬ 
dictions for phase-resolved spectra which can be tested 
by observations. It has been applied to four magnetars 
with available phase-resolved data above 10 keV, and in 
all cases t he model was found consistent with the data 
(lAn et al.l 120131 |Hasc oet. Beloborodov fe den Hartogl 

12014 IVogel ei~aLll2014ll. ~ 

The magnetar IE 1841—045 has a surface dipolar 
magnetic-field strength of B = 3.2 x 10 19 (PP) 1 / 2 G = 
6.9 x 10 14 G, estimated from the spin period and the 
spin-down rate of P = 11.8 s and P = 4 x 10~ n s s _1 , 
respectively, assuming the standard vacuum dipole spin- 
down model. The source has previously shown occasional 
X-ray bursts with ener gi es of ~ 10 38 erg a s suming a 
distance of 8.5 kpc (IKumar fe Safi-Harbll2010t ILin et all 

120 1 llfDib fe K as nil 201 41) . No tails or significant enhance¬ 

ment in the persistent emission were ob served following 
the bursts. IKumar fe Safi-Harbl (120101) reported detec¬ 
tion of emission lines >20 keV in the burst spectrum mea¬ 
sured with the Swift Burst Alert Telescope (BAT). Fur¬ 
thermore, the authors argued that the source brightened, 
and the emission b e came softer after the burst activity. 
However. ILin et ahl (1201 lh did not find any emission line, 
flux enhancement or spectral change s at the burst epoch 
i n con tradiction with the results of IKumar fe Safi-Harbl 

dma. 

IE 1841—045 is one of the brightest magnetars in the 
hard X-ray band ab ove 10 keV. Its pers istent emission 
has b een studied by iKuiper et al.l (|2006T) and lAn et all 
(1201311 . The measured photon index in the hard X-ray 
band is T ~ 1.3, and the puls ed fraction was reported 
to increase with photon energ y. lAn et al.l (12 0 1311 ap plied 
the coronal outflow model of iBeloborodovI (T2013D to a 
50-ks NuSTAR observation and constrained the hard X- 
ray emission geometry to two possible solutions. Fur¬ 
thermore, an interesting feature was found in the pulse 
profile in the 24-35 keV band, which may be associated 
with a spectral feature in this band. 

In this paper, we further investigate the properties of 
persistent emission of IE 1841—045 using a new 350-ks 
NuSTAR observation and 3-year monitoring observations 
by Swift. Fortuitously, the source was actively bursting 
during the NuSTAR observation, which provided an op¬ 
portunity to study the bursts in addition to the persistent 
emission. We describe the NuSTAR and archival Chan¬ 
dra ., XMM-Newtor and Swift observations used in this 
paper in Section [2] and present the results of the NuS¬ 
TAR data analysis in Section [3] We then present the 
Swift monitoring observations and the data analysis in 
Sections 0] and 14.11 We discuss the results of data analy¬ 
ses in Section [5] and summarize our main conclusions in 
Section [6] 


IE 1841—045 was observed with NuSTAR 
(jHarrison et al.l I2013H between 2013 September 5 
and 21 in a series of exposures with durations 40-100 
ks with a total net exposure of 350 ks. Two X-ray 
bursts were detect ed with the Fermi Gamma-Ray Burst 
Monitor (GBM) dCollazzi. Xiong. &; Kouveliotoul 12014 

iPal’shin et al.ll2014fl on 2013 September 13. during the 

NuSTAR observation period. Fortunately, the bursts 
were also recorded in the NuSTAR data. In addition 
to the Fermi reported bursts, NuSTAR serendipitously 
detected s ever al more bursts. We report on the bursts 
in Section I3TT1 To study the persistent emission of the 
source, we also analyzed archival observations made 
previously with NuSTAR, Chandra, XMM-Newton and 
Swift to have better statistics and to constrain the 
persistent properties of the source in the soft band. The 
observations used in this study are listed in Table [Q 
Note that all the soft-band observations and the first 
NuSTAR observatio n (Obs. ID 30001025002) were 
reported previously (lAn et al.ll2013l A13 hereafter). 

The NuSTAR data were processed with nupipeline 
1.3.1 along with CALDB version 20131223. We used 
default filters except for PSDCAL for which we used 
PSDCAL=N0. The PSDCAL filter is for laser metrology cali¬ 
bration of the relative positions of the optics and detec¬ 
tors. This observation was affected by times when the 
metrology laser went out of range. By not using the de¬ 
fault, the pointing accuracy may be slightly degraded, 
but good time intervals (GTI), and hence exposure time, 
increaseQ We note this situation is unusual and specific 
to this observation. We verified that the analysis results 
described below do not change depending on the PSDCAL 
filter setting. However, we note that some of the burst 
data could be retrieved only with PSDCAL=N0. 

We also analyzed archival XMM-Newton and Chandra 
observations (see Table[T|). The XMM-Newton data were 
processed with Science Analysis System (SAS) 12.0.1, 
and the Chandra data were reprocessed using chan- 


A See http://heasarc.gsfc.nasa.gov/docs/nustar/analysis/mistar_swguide.pdf 
for more details 
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dra_repro of CIAO 4.4 along with CALDB 4.5.3. We fur¬ 
ther processed the cleaned data for analysis as described 
below. Uncertainties below are at the lcr confidence level 
unless stated otherwise. 

3. DATA ANALYSIS AND RESULTS FOR THE NUSTAR 
OBSERVATIONS 

In this Section, we present data analysis results for 
bursts and persistent emission measured with the obser¬ 
vations in Tabic [T] (Sections I3.H3.3D . We fit the persis¬ 
tent and phase-resolved spectra with the coronal outflow 
model and show the results in Section [XU 

3.1. Burst Analysis 

3.1.1. Temporal and Spectral Properties of the Bursts 

We performed a comprehensive search for bursts in 
the NuSTAR light curves. We extracted event time 
series, applied the barycenter correction using the po¬ 
sition R .A.=18 h 41 m 19 s .343, Decl.=-04°56'11"8 (J2000; 
iDurantl 120051) . binned the light curves with a variety of 
bin sizes ranging from 0.01 to 10 s, and searched for 
time bins which contained more counts than expected 
above background including source persistent emission 
using Poisson statistics. The background was extracted 
in an interval 10-pulse periods long using the same ex¬ 
traction region, just prior to the time bin being con¬ 
sidered. In total, we found 7 time bins which are sig¬ 
nificantly above the mean level after considering the 
number of trials. Note that two of the seven signif¬ 
icant bins turned out to be produced by one burst 
(burst 5; see below), hence we found six bursts dur¬ 
ing our observations. The significance of the bursts is 
high (p < 10~ 10 ) but only on short timescales, e.g., 
^0.1 s. We list the burst times in Tabic [21 Note that 
bursts 5 and 6 were p reviousl y reported based on Fermi 
GBM detections dCollazzi. Xiong. &; Kouveliotoul 12014 

iPal’shin et al.ll201fl). 

We note that some of the bursts may not be fully sam¬ 
pled due to high count rates. Specifically, the maximum 
count rate of NuSTAR detectors is ~400 cps limited by a 
deadtime of 2.5 ms per event. In addition, an event that 
was detected in a very short elapsed time from the previ¬ 
ous event is regarded as background and filtered out dur¬ 
ing the standard pipeline process. We investigated these 
effects by looking into the elapsed time for each event, 
and found that those for events in two high-count time 
bins were significantly shorter than in other time inter¬ 
vals. We further reprocessed th e observation data w ith 
a relaxed elapsed-time filter (see lMadsen et al.ll2015l for 
more details) and were able to recover an additional 58 
events in a R = 120" circular aperture in the 3-79 keV 
band for the two time bins combined. From this study, 
we find that the two high-count bins were actually con¬ 
nected in time (i.e., the gap between the two bins was 
filled by the recovered events) and the livetime of the 
detectors was ~l/300 of the exposure. We also investi¬ 
gated other high-count time bins using a relaxed filter, 
and were able to recover marginally additional events 
only for burst 6 (10±8 events), the other GBM detected 
burst. Below, results for burst properties are obtained 
with the data processed with the relaxed filter. 

Although it was reported that the Fermi-detected 
bursts are likely from the magnetar IE 1841—045 


(iPal’shin ct al.l[2014f) . the localization was not unambigu¬ 
ous. In order to localize the bursts better and see if they 
are really produced by the magnetar IE 1841—045, we 
used the NuSTAR data to measure the position of bursts 
2, 3, 5, and 6, which had sufficiently many counts for such 
an analysis. We projected their 3-79 keV event distri¬ 
butions collected for 2 s onto R.A. and decl. axes, and 
fit the projected profiles with a Gaussian plus constant 
function. The results for the burst location offsets from 
the IE 1841-045 position were A R = 9" ± 2", 8" ± 2", 
9" ± 3", and 5" ± 2" for bursts 2, 3, 5 and 6, respec¬ 
tively, where the quoted uncertainties are purely statis¬ 
tical (lcr). Note that aspect reconstruction accuracy of 
NuSTAR is 8" (90%), and so the measured positions are 
consistent with that of the magnetar IE 1841—045. 

In order to characterize the properties of the bursts, we 
fit the short-term (~10-s) light curves around the bursts 
with exponentially rising and falling functions 


m 


_ f Ae { - t - To '>/ T ' + Ci i<T 0 , 

~ \(A-G 2 )e-( t - T “)/ Tf +Ci + C 2 t > To, 


( 1 ) 


where A is the amplitude, To is the burst peak time, T r 
is the rising ti me, Tf is the fa lling time, and C\^ are 
constants (see iGavriil, Dib. fe Kasnil I2011L for a differ¬ 
ent model). Since there are only 10-120 counts detected 
within a 1-s window around each burst, we extracted 
events from the whole detector and used a maximum 
likelihood optimization without binning. The peak count 
rates are very high and the T 90 durations are ^1 sec. We 
present the results in Table [2] and show the bursts mor¬ 
phologies in Figure [T] Some bursts occurred within 1 sec 
of each other. Note that we do not find a clear increase 
in the tail emission except for in burst 5. 

The spectrum of a burst can provide information on the 
burst mechanism. Therefore, we extracted ~0.2-s spec¬ 
tra in circles with radius R = 120" centered at the source 
position for bursts 2 and 5, and fit the spectra with sin¬ 
gle component models; we tried both a blackbody and 
a power-law model. In order to remove the persistent 
emission, we extracted background counts in the same 
region as we used for the source but in a 1-ks pre-burst 
interval. We did not attempt to fit the burst data with a 
multi-component model because there were too few pho¬ 
tons collected during the 0.2-s intervals. We grouped the 
spectrum to have 1 count per spectral bin, and used Istat 
dLorcdo 1119921) in XSPEC 12.8.lg. The NuSTAR bandpass 
is not sensitive to the relatively low Ah of the source, and 
we set it to 2.05 x 10 22 cm -2 which we obtained by joint¬ 
fitting of the soft-band spectra with a broken power-law 
model (see Section 13. 2. 3D . We use this value and the 
tbabs absorption model in XSPEC throughout this pa¬ 
per unless noted otherwise. The burst spectra can be 
described with a power-law model having T = 1 — 2 or a 
blackbody model with kT = 3 — 5 keV. The results of the 
fits are shown in Tabled Note that we did not detect the 
hig h temperature blackbody co mponent (kT =13 keV 
ICollazzi. Xiong, fc Kouvcliotoull2014f) , which is probably 
because of the low statistics at high energy. 


3.1.2. Burst Tails 

Short X-ray bursts from mag netars can exhi bit long 
emission tails, lasting for hours (lAn et al.l I2Q14T) . In or¬ 
der to search for long tails, we binned the light curves into 
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Table 2 

Deadtime-corrected Properties of TVuS'TA.R-detected Bursts from IE 1841—045 


Burst 

To 

0“ 

T r 

Tt 

A 

Ci 

c 2 

^90 b 

-^evt ^post —pre 


(day) 


(s) 

(s) 

(cps) 

(cps) 

(cps) 

(s) 

(cts) (cts) 


1 

0.35836545 

0 2266~*~ U ' UUU3 

U.ZZDD_ 0 0002 

0 007 +u ' uul 

u.uu/ _ 0 001 

0 047 +u ' ul7 

— Q. Q i 9 

17401 '™ 

12ta 

< 16 

0 i2 +U (J4 

U - 1Z -0.05 

5 


2 

0.35837490 

0 2958" l "°' c ^ 2 

U.^300_g QQQ2 

0 006" l ’ < ^ li2 
u.uuo_ 0 002 

0 09 +0 ' 01 
u.uu_ 0 01 

700t “° 

5 . 511 “ 

5 9+ 1 - 8 

ri oq + 0.03 

U.ZZ _ 0 03 

61 

75 ± 41 

3 

0.60981692 

0 5 58 0002 

U.OOOO_Q 0001 

0 . 014 ^ 0.001 

0 . 011^:002 

2010j :f° 0 

121 ' 

c+ 3 

°-3 

U.UO ( _o .007 

27 

16 ± 40 

4 

7.27821493 

0 4463 +0 ' 0002 

U.‘±‘±UO_Q >0002 

o.oosslD™™ 

0 . 052^:02 

840 ± 1 ™ 

161 “ 

< 3.4 

f) 1 O + 0.05 
U.io_ 0 Q5 

11 

26 ± 40 

5 

8.62801288 

0 1 075 - * -0 ' 0001 
u.iu/o_ 0 0001 

0 0090’ 1 ’ 0 0003 
u.uuuu_ 0 000 3 

n m q^+o.0004 
O.Ol84_ 0 0004 

67000 

isi; 

< 3 

0 0631+0 000 7 

U.UOOl_ 0 QQQ7 

22 

211 ± 43 

6 

8.759684723 

n ooc'7Q - t _ 0 -00001 
U.ocJb ( ^- 0.00001 

< 0.0006 

0 0249"* -0 ' 0016 

u.uz4iy_ 0 0018 

80001 )™“ 

14l ) 

< 4 

< 0.059 

31 

7 ±41 


Notes. Parameters for the short timescale light curve. To is the burst arrival time and is days since MJD 56540 (barycentric dynamical 
time). T r f are the rising and falling times for the burst light curves, A is the peak count rate, C 1.2 are the constant level of the light curves 
before and after the bursts (see Equation [TJ, T 90 is the time interval which includes 90% of the burst counts estimated with the exponential 
functions, N ev t is the number of events within T 90 , and A pos t— pre is the difference in numbers of photons contained in the pre- and the 
post-burst 200 -s intervals. 

a Spin phase corresponding to Tp, w here phase zero is defined at the pulse minimum (T re f = 56540.32899020 MJD), the same as that for 
the timing analysis in Section l3.2.2l 

b Since T 90 for the whole burst is not always well defined because the constants Ci ,2 are different before and after the burst peak, Tgo’s 
for the rising and the falling function were calculated separately and then summed to obtain that for the burst. When only an upper limit 
is available for T r or Tf, we used the upper limit to calculate T 90 and show it without uncertainties. 
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Figure 1. Deadtime-corrected light curves of the bursts in the 3—79 keV band. 


Table 3 

Best-fit parameters of the burst spectra for a ~0.2 s interval 
around the burst peak 


Burst 

r 

3-79 flux 8 

10 -8 erg s -1 cm -2 

Istat/doi 

2 

1.6(3) 

3.5(1.2) 

37/46 

5 

0.96(24) 

1300(600) 

53/56 

Burst 

kT 

(keV) 

E 

^BB 

10 38 erg s 1 

Istat/doi 

2 

3.3(3) 

1.5(3) 

48/46 

5 

4.8(5) 

400(100) 

59/56 


a Deadtime-corrected flux. 

b Deadtime-corrected bolometric luminosity for an assumed dis¬ 
tance of 8.5 kpc. 

20-s bins, and compared counts in 10 bins before and af¬ 
ter a burst, excluding the burst bin. We then compared 


the pre- and post-burst counts, and found that the differ¬ 
ence is significant only for burst 5 (A pos t_p re = 211 ± 43 
counts for a 200-s time interval). We show the light 
curves for the bursts in Figure[2l and report the results 
in Table [2] We performed the same study on different 
timescales (e.g., 2 s and 50 s), and found the same results; 
the difference is significant only for burst 5 
In order to search for spectral evolution after burst 5, 
we extracted spectra within a radius R = 120" in the tail 
of burst 5 in several time intervals excluding the burst 
(T > T 0 + 0.5 s; see Table O. Each interval had > 100 
events above the persistent emission plus background. 
We then fit the spectra with a blackbody and a power- 
law model. We show the results of the power-law fit in 
FigureOH The spectral shape did not change significantly 
over ~2-ks of tail emission, and the 3-79 keV flux decay 
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Figure 2. Observed long-term light curves of the bursts in the 3-79 keV band. Note that bursts 1 and 2 are very close in time (^5l s), 
and that the light curve for burst 4 is not shown because the burst was not detected in a 20-s time bin. The blue dotted line shows the 
pre-burst emission level, and the red line shows the scaled Good Time Interval (GTI). 


is well described with a power-law function, having a 
decay index of 0.45 ± 0.10. Note that the decay index 
was measured after taking into account the covariance 
between the flux and the spectral index in the fit. A 
similar decay index is measured for the lumin osit y when 
using a blackbody model. Note that lAn et all (12014H 
measured the flux evolution index in the tail to be 0.8- 
0.9 for IE 1048.1—5937, much steeper than our measure¬ 
ment for IE 1841—045. The count enhancement at later 
times seems to be spiky, and might have been caused 
by undetected activity there. Since the other bursts do 
not have significant tail emission, we were not able to 
measure their spectral evolution. 

The tail spectrum of burst 5 is soft compared to the 
burst spectrum, and may be similar to the persistent 
emission. In order to see if the spectra of the burst 
tails are significantly different from the persistent source 
emission, we studied the persistent emission in a pre¬ 
burst interval in which approximately the same number 
of events was collected as in the tail spectra. We ex¬ 
tracted a persistent spectrum in a 100-s pre-burst time 


interval, modeled it with a single component model, ei¬ 
ther a power law or a blackbody. The persistent spec¬ 
trum over the short interval was well described with sin¬ 
gle component models, and we find the best-fit parame¬ 
ters are T = 2.31 ± 0.24 or kT = 1.58 ± 0.13 keV, similar 
to the tail spectrum. We show the persistent levels with 
blue dashed lines in Figure 4. 

Since the spectral shape did not change significantly 
over the tail interval, we fit the combined (~1.8 ks) tail 
spectrum to a power-law or a blackbody model after re¬ 
moving the persistent emission. The spectrum is well fit 
by a power-law model (y 2 /dof=122/146) with a photon 
index T = 2.2±0.2 (green lines in Figurc[3]left), similar to 
the 100-s persistent spectrum. A blackbody model also 
fits the data with kT = 2.1±0.2 keV (y 2 /dof=137/146). 
We tried to fit the combined tail emission with the model 
for the persistent spectra obtained below (Section 13.2.31) 
using the same fit parameters except for the normaliza¬ 
tion constants (Table [4|. The model fits the spectrum 
well (x 2 /dof=171/147 with the null hypothesis probabil¬ 
ity p = 0.09). However, we see a trend in the fit residuals 
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Figure 3. Evolution of the persistent-emission-removed spectrum of the tail of burst 5 measured using a power-law model (left) and the 
integrated spectra (right). We used events collected after T = To+0.5 s to remove the burst emission. Flux is in units of 10“ 11 erg s —1 cm -2 
in the 3-79 keV band. The red dotted line in the left panel shows the power-law function that best describes the flux evolution, blue dashed 
lines show the lcr range of the quantities measured in a 100-s interval before the burst. The green dashed lines in the left bottom show 
the lcr range of the best-fit power-law index for the integrated tail spectrum. The integrated tail spectrum with a blackbody plus broken 
power-law fit with spectral shape parameters frozen at the values in Table [4] is shown in the right panel. 
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which suggests that the tail spectra may be slightly softer 
than the persistent one (see Figure [3] right). 

3.2. Persistent Emission 

In order to study the persistent emission, we removed 
the burst intervals using time filters. We used 20-s win¬ 
dows centered at the burst peak times for all the bursts 
except for burst 5 for which we used a 2-ks window 
because of the tail. For the soft-band spectrum below 
10 keV, we used the Chandra (only Obs. ID 730 due 
to pile-up in Obs. ID 6732), XMM-Newton , and Swift 
data (see Table [TJ, the same observations as used by 
A13. Although these observations were taken long be¬ 
fore the NuSTAR observation, we show below that the 
source emission properties have been stable over 10 years 
(see also A13). For the NuSTAR data, we used a circu¬ 
lar aperture with R = 60" for the source and an annu¬ 
lar aperture with inner and outer radii of R = 60" and 
R = 100" for the background, respectively. Note that 
~15% of the source events fall in the backgro und region 
because of the NuSTAR PSF i An et al.ll2014H . However, 
not all the ~15% of the source flux is subtracted as back¬ 
ground in the spectral fitting since we scale the area of 
the background region to that of the source for spec¬ 
tral fits, and ~10% of source events will be lost during 
background subtraction. We take into account this effect 
using a normalization constant. 

3.2.1. Timing Analysis 

For our timing analysis, we extracted source events 
from the new observations in the 3-79 keV band within 
a radius R = 60" of the nominal source position, and 
divided the events into subintervals consisting of ~5000 
counts each. Note that we have used different subin¬ 
tervals and found that the results do not change signifi¬ 
cantly. Each subinterval was then folded at the nominal 
pulse period to yield pulse profiles each with 64 phase 
bins. In order to perform phase-coherent timing, we 
first cross-correlated the pulse profiles to measure the 
phase for each subinterval. We then fit the phases to 
a quadratic function using frequency (v) and its first 
derivative (z>), <f{t) = (f 0 + v(t — T re f) + v(t — T re f) 2 /2, to 
derive a timing solution. We produced a high signal-to- 
noise ratio template by coherently combining the pulse 
profiles using the timing solution, and cross-correlated 
the pulse profiles with the template in order to refine 
the timing solution. We show the residuals after the fit 
in Figure [I] We find that the best timing solution dur¬ 
ing the observations has parameters P = 11.79234(1) s 
and P = 4.2(2) x 10 -11 s s -1 , implying a magnetic field 
strength of 7 x 10 14 G. Note that we did not include the 
first NuSTAR observation (Obs. ID 30001025002) in this 
study because of phase ambiguity. We verified the timing 
solution by measuring the period for the individual ob¬ 
servations including the first NuS TAR observation (Obs . 
ID 30001025002) using the U-test (|de Jaeer et al.lll989lf . 
and by fitting the measured period to a linear function of 
period evolution P{t) = Pq + P(t — T re f) s, which yields 
P = 11.792344(4) s and P = 4.05(9) x 1CT 11 s s" 1 . We 
find that the results of the two methods are consistent 
with each other and with the results of our Swift moni¬ 
toring (see Section |4~T1) . 


i.or 



-I.QLl 


0 5 10 15 20 

Time since MJD 56540.3 (Days) 

Figure 4. Timing residuals after fitting the pulse phases in the 
3-79 keV band for observations 30001025004-12. 

3.2.2. Pulse Profiles and Pulsed Fraction 

The pulse p rofile of IE 1841—0 4 5 is known to c hange 
with energy (IKuiper et all 120061 lAn et all I2013T) . In 
particular, A13 found that the pulse shape in the 24- 
35 keV band is different from those in the adjacent en¬ 
ergy ranges, which suggested the existence of a spectral 
feature. However, no firm conclusion could be made due 
to limited statistics. We investigate this here with much 
better statistics. 

We produced pulse profiles for individual observations, 
and aligned them with the template. Backgrounds were 
extracted from an annular region with inner and outer 
radii of 60" and 100" and were subtracted from the pulse 
profile of the source region. We verified that the pulse 
shape has not changed significantly over the ~300 days of 
NuSTAR observations (Table [TJ by comparing the pulse 
profile of individual observations with the combined pro¬ 
file in several energy bands (e.g., 3-6 keV, 6-10 keV, 10- 
15 keV, and 15-79 keV). We show the combined pulse 
profiles in several energy bands in Figure [S] Note that 
we find double-peaked structure similar to that seen by 
A13, e.g., in the 17 -33 keV band (see Figurc[S]). However, 
with the much better statistics we have now, we find that 
the pulse shape does not change suddenly but instead 
changes gradually with energy. This does not support 
the existence of a narrow spectral feature as suggested 
by A13. 

We find that the pulse shape at higher energies be¬ 
comes more complicated, sometimes showing three peaks 
(e.g., 33-38 keV in Figure [5] a new peak seems to ap¬ 
pear at phase ~0.5). In order to see if the triple-peaked 
structure at high energies (^ 25 keV) is significant, we 
fit each pulse profile with a harmonic function in which 
the number of harmonics contained was varied between 
zero (constant) and three. In the fit, we calculate the y 2 
value and the F-test probability by adding higher-order 
harmonic functions one by one. From this study, we find 
that the pulse profiles below 38 keV are generally well fit 
with sum of two harmonics, and the others with a single 
harmonic function; adding one more harmonic to these is 
not statistically required (99% confidence). We further 
fit the pulse profiles with a sum of the first two har¬ 
monics plus a fifth harmonic because the triple-peaked 





























































structure is best described with a fifth harmonic. The 
F-test probability for adding the fifth harmonic shows 
an improvement to the fit with 99.7% confidence in the 
45-55 keV band and with 98.6% confidence in the 33- 
38 keV band. However, these may not imply a signifi¬ 
cant detection when considering the number of double- 
peaked profiles we have. Therefore, we conclude that the 
double-peaked structures are statistically significant but 
the triple-peaked ones are only marginally so. 

The pulsed fraction of the source has pr evious l y been 
meas ured to be increasing with energy (jKuiper et al.l 
|200& A13). Note that A13 did not attempt to measure 
the area pulsed fraction defined by 


PF — 

1 x area 


Ei^* Pmin) 

E iPi 


( 2 ) 


where pi is counts in i-th phase bin and p m in is the counts 
from the lowest bin in the pulse profile, due to insufficient 
statistics, and that PF area is known to be biased upwards 
in the low counts regime. Since we now have much bet¬ 
ter statistics, we measured the area pulsed fractions in 
different energy bands and show them in Figure [5] Even 
with good statistics, the area pulsed fraction is known 
to be biased upwards (see Appendix), and so we show 
two alternative measures of the pulsed fraction as well. 
The first alternative is to fit the pulse profile with a har¬ 
monic function, and use the best-fit function to calculate 
the pulsed fraction (PEgt), which will remove the bias 
caused by incorrect identification of the baseline level by 
selecting the minimum phase bin. We find that the pulse 
profiles are well described with two harmonics having 
X 2 /dofi"4 (Figure [SJ except for the lowest energy band 
for which the two-harmonic function yields an unaccept¬ 
able fit, with p = 10 -4 . The large y 2 in this case is 
mostly due to a sharp step in the pulse profile at phase 
~0.2, but the overall pulse shape agrees with the best- 
fit function well. Furthermore, rebinning the phase can 
make the fit acceptable. Therefore, we used two har¬ 
monics for the fit function, calculated PF area using the 
best-fit parameters, and show the pulsed fractions in Fig¬ 
ure [G] (blue squares). We also calculated the rms pulse 
amplitude defined by 
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similar to that seen by A13 (their Figure 3) if we use dif¬ 
ferent energy bins. However, the overall trend is similar; 
we do not see any rapid increase in the pulsed fraction 
with ene rgy above lO keV. This is different from previous 
reports (iKuiper et aJl2()()( >; of pulsed fraction increasing 
with energy and approaching 100% at 100 keV. 

3.2.3. Phase-Averaged Spectral Analysis 

Since it is possible that the source has different emis¬ 
sion properties during the bursting periods, we compared 
the spectral properties of individual observations. In this 
study, we did not include the soft-band spectra because 
they were taken at much earlier epochs. We jointly fit 
all the NuSTAR spectra in the 6-79 keV band with an 
absorbed broken power-law model in order to minimize 
the effect of the blackbody component, which is negligi¬ 
ble above 6keV, and found that the spectral shapes for 
the six NuSTAR observations (see Table[T]) are consistent 
with one another. 

Since there is no clear evidence that the source spectral 
shape has varied during the NuSTAR observations, we 
used all the NuSTAR observations for the phase-averaged 
spectral analysis. Furthermore, we used the soft-band 
data as well in this study since the shap e of the soft- 
band spectrum is also kn own to be stable dZhu fe Kaspil 
120101 : iDib fc Kaspil l2014lh We tied all the model pa¬ 
rameters between NuSTAR , Swift , XMM-Newton, and 
Chandra except for the cross-normalization factors. The 
normalization constant for NuSTAR FPMA (Obs. ID 
30001025002) was set to be 0.9 as a reference in order to 
account for the source contamination in the background 
region (Section 13.2.21) . To fit the data, we grouped the 
spectra to have at least 20 counts per bin. We used 
an absorbed blackbody plus broken power law, and an 
absorbed blackbody plus a double power law to fit the 
3-79 keV NuSTAR data and the 0.5-10 keV soft-band 
data. We present the results in Table []] and the spectra 
in Figure [7] 

We note that the spectral parameters we report in Ta¬ 
bleware slightly different from those reported previously 
(A13). In order to see if the difference is due to the up¬ 
dated calibration, we analyzed the same data set that 
A13 used (Obs. ID 30001025002), and were able to re¬ 
produce their results except for the flux. The flux we 
measure is lower by ~15% than what A13 reported, be¬ 
cause of a ~15% increase in the NuSTAR effective area 
from CALDB 20131007Q The hard-band component is 
much better constrained with the new long exposures, 
and thus the new results we report are more accurate. 
We note that the new parameters in Table [4] are not in¬ 
consistent with the data A13 used, providing acceptable 
fits to the data with x 2 /dof =2965/2878 and 2914/2878 
for the blackbody plus broken power law and the black¬ 
body plus double power law, respectively. 


is the Fourier power produced by the noise in the data, pt 
is the number of counts in ith bin, N is the total number 
of bins, (7p i is uncertainty in pi, and n is the number 
of Fourier harmonics included, in this case, n = 6 (see 
lArchibald et al.112014 for more details). We show the 
results in Figure [6] (red diamonds). 

We note that the small-scale features in Figure [6] 
change when we use different energy resolution. For ex¬ 
ample, a sudden jump sometimes appears at ~30 keV 


3.3. Phase-Resolved and Pulsed Spectral Analyses 

We conducted a phase-resolved spectral analysis for 
ten phase intervals to study distinct features in the pulse 
profiles (see Figure [5] for pulse profiles). We did not use 
the Swift XRT or XMM-Newton MOS data in this study 
because their temporal resolutions were insufficient. The 

J http://heasarc.gsfc.nasa.gov/docs/heasarc/caldb/nustar/docs 
/release_20131007.txt 
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Figure 5. Background-subtracted pulse profiles for IE 1841—045 measured with NuSTAR in various energy bands. The average value is 
shown in a black dotted line and a Fourier reconstructed profile with five harmonics is shown in red in each panel. 


Table 4 

Phenomenological spectral fit results for IE 1841—045 


Phase 

Data 3- 

Energy 

(keV) 

Model 3 

JV H 

(10 22 cm' 2 ) 

kT 

(keV) 

F c 

4 S 

Fbreak/Fs d 

(keV/ ) 

r h //3 e 

Fh* 

Lbb s 

X'Vdof 

o 

1 

o 

6 

N,S,X,C 

0.5-79 

BB+BP 

2.05(3) 

0.491(5) 

1.95(1) 

13.5(2)/-•■ 

1.24(1) 

5.88(6) 

1.64(5) 

8060/7930 

o 

l 

o 

d 

N,S,X,C 

0.5-79 

BB+2PL 

2.49(5) 

0.443(9) 

2.82(8) 

■ - - /1.53(6) 

0.97(4) 

4.70(6) 

1.15(9) 

7931/7930 

Pulsed 

N,X,C 

0.5-79 

PL 

2.05 h 




1.83(3) 

1.4(1) 


1236/1749 

Pulsed 

N 

3-79 

PL 

2.05 




1.81(3) 

1.4(1) 


1128/1591 

Pulsed 

N 

5-79 

PL 

2.05 




1.73(4) 

1.4(1) 


806/1140 

Pulsed 

N 

10-79 

PL 

2.05 




1.47(7) 

1.5(2) 


349/510 

Pulsed 

N 

15-79 

PL 

2.05 




1.28(12) 

1.6(3) 


174/279 

Pulsed 

N 

20-79 

PL 

2.05 




1.0(2) 

1.5(4) 


113/159 


a N: NuSTAR, S: Swift, X: XMM-Newton, C: Chandra. 

b BB: Blackbody, PL: Power law, BP: Broken power law, and 2PL: Two power laws in XSPEC. 
c Photon index for the soft power-law component. 

d Break energy for the BB+BP fit or soft power-law flux in units of 10 -11 erg s -1 cm -2 in the 3—79 keV for the BB+2PL fit. 
e Photon index for the hard power-law component. 

1 Flux in units of 10 —11 erg cm -2 s —1 . The values are only the power-law (hard power-law) flux in the 3—79 keV band for the BP (PL, 
2PL) model. 

g Blackbody luminosity in units of 10 35 erg s —1 for an assumed distance of 8.5 kpc IITian fe f.calivll2008ii . 
h A+ for the pulsed spectral analysis was frozen. 
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Energy (keV) 

Figure 6. Pulsed fractions at several energy bands measured with 
NuSTAR. Black triangles : area pulsed fraction (Equation [2j; blue 
squares : area pulsed fraction measured using harmonic fit; and red 
diamonds : rms pulsed fraction (Equation [3}. 
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Figure 7. Phase averaged spectra of IE 1841—045 and the fit re¬ 
sult. Chandra , XMM-Newton, and Swift data cover below 10 keV, 
and NuSTAR data cover 3-79 keV (see Table [J| for observation 
summary). Each component of the best-fit model, an absorbed 
blackbody plus double power law, is shown in lines. See Table [4] 
for best-fit model parameters. 

Chandra and XMM-Newton PN data were phase-aligned 
with the NuSTAR data by correlating the pulse profiles. 

We binned the NuSTAR and soft-band instrument 
spectra to have at least 20 counts per spectral bin, and 
froze the cross-normalization factors to those obtained 
with the phase-averaged spectral fit. We fit the spectra 
with the two models that we used for fitting the phase- 
averaged spectrum: an absorbed blackbody plus broken 
power law and an absorbed blackbody plus double power- 
law model. We find that both models explain the data 
well, having y 2 /doKl.003 for dofs of 1413-1966. The 
spectra vary with spin phase, having harder power-law 
spectra when the flux is high. However, the detailed vari¬ 
ation depends on the spectral model used. We show the 
results in Figure [8] 

In order to see if there is a spectral feature that 
shifts with e nergy as was seen in SGR 0418+5729 
(|Tiengo et alJ 12013 1. we produced an energy-phase im¬ 
age using 25 phase bins and 40 energy bins in the 3- 
79 keV band. We first divided counts in each pixel in the 
energy-phase 2-D map by the phase-integrated counts in 
the same energy bin and present it in the left panel of 


Figure 0 which shows similar structures to the energy- 
resolved pulse profiles (Figure [5]). We then divided the 
map further by the energy-integrated counts in the same 
phase bin (Figure 0 in order to have a better contrast, 
and find no clear phase-dependent feature in the image 
(Figure 0 right). We also tried different binning and 
found the same results. 

We measured the pulsed spectrum in the 0.5-79 keV 
band in order to see if it is significantly harder than 
the phase-averaged sp ectrum as see n in other hard X- 
ray bright magnetars dKuipcr et al.1 120061) . We grouped 
the spectra to have at least 200 counts per spectral bin 
and subtracted the spectrum in the phase interval 0.9-1.1 
(the DC l evel) f rom the phase-averaged spectra obtained 
in Section l3.2.3l We then jointly fit the broad-band spec¬ 
trum with a power-law model letting the cross normal¬ 
ization constants vary. We used 1st at and y 2 statistics 
and found that they give consistent results. 

A simple power-law model with a photon index of 1.8 
fits the 0.5-79 keV data well (reduced x 2 /dof<l). In or¬ 
der to verify the measurements for the pulsed spec trum 
above 15 keV (T = 0.72 ± 0.15. iKuiper et ahll2006lf . we 
restricted the fit range to high energies. As the lower en¬ 
ergy cutoff is increased, the power-law index decreases, 
consistent with spectral hardening. The results are sum¬ 
marized in Table 0 

We also estimated pulsed fractions in the hard band us¬ 
ing the spectra (defined as the ratio of pulsed and total 
flux densities) in order to compare with those in Figure0 
We fit the total and the pulsed spectra (^15 keV) to sin¬ 
gle power-law models. The power-law index of the total 
spectrum above 15 keV is 1.19 ± 0.02 (1.13 ± 0.02 above 
20 keV), slightly smaller than what we obtained using the 
absorbed blackbody plus broken power-law model (see 
Table 0. The power-law index of the pulsed spectrum 
above 15 keV is 1.28 + 0.12 (1.0 + 0.2 above 20 keV), sim¬ 
ilar to that of the total spectrum. This suggests that the 
pulsed flux does not rapidly increase in the hard band, 
as also seen in Figure 0 

3.4. Spectral fits with the e^ outflow model 

A13 found that the properties of the persistent hard 
X-ray emission of IE 1841—045 wa s consistent with th e 
coronal outflow model proposed bv iBeloborodovl (1201311 . 
The model envisions an outflow of relativistic electron- 
positron pairs created by electric discharge near the neu¬ 
tron star. The outflow fills the active “j-bundle” — a 
bundle of closed magnetospheric fi eld lines that carry 
electric current (IBeloborodovl 1200911 . The pair plasma 
flows out along the magnetic field lines and gradually 
decelerates as it scatters the thermal X-rays. It radi¬ 
ates most of its kinetic energy in hard X-rays before the 
e ± pairs reach the top of the magnetic loop and anni¬ 
hilate. The magnetic dipole moment of IE 1841—045 is 
estimated from its spin-down rate, /i« 7 x 10 32 G cm 3 , 
assuming t he neutron star radius to be 10 km. Simila r 
to A13 and lHascoet, Beloborodov fe den Hartog I (I2014T) . 
we assume a simple geometry where the j -bundle is ax- 
isymmetric around the magnetic dipole axis. However, 
instead of assuming that the j-bundle emerges from a 
polar cap, its footprint is allowe d to have a r ing sh ape. 
This possibility was introduced in lVogel et al.l (I2014H . be¬ 
cause the NuSTAR data for IE 2259+586 favored a ring 
footprint over a polar cap. The model has the following 
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Figure 8. Results of our phase-resolved spectral analysis. Blackbody luminosity (Lbb) is i n units of 10 35 erg s 1 . and power-law flux 
(F pl ) is in units of 10“ 11 erg s~ 1 crrO 2 in the 3-79 keV band. Gamma), is the power-law index of the hard power-law component and 
Gamma s is for the soft component. 
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Figure 9. Background-subtracted energy-phase count images in the 3-79 keV band produced using 25 phase bins and 40 energy bins. 
Counts in each pixel were divided by the phase-integrated counts in the same energy bin (left), and then by the energy-integrated counts 
in the same phase bin (right). Two phases are displayed for clarity. 


parameters: (1) the power Lj of the e ± outflow along 
the j-bundle, (2) the angle a mag between the rotation 
axis and the magnetic axis, (3) the angle /3 0 bs between 
the rotation axis and the observer’s line of sight, (4) the 
angular position 9j of the (/-bundle footprint, and (5) the 
angular width A 6j of the j-bundle footprint. In addition, 
the reference point of the rotational phase, </>o, is a free 
parameter, since we fit the phase-resolved spectra. 

We foll ow the method pre sented in 
lHascoet. Beloborodov fe den Hartog I (12014D and explore 
the whole parameter space by fitting the phase-averaged 
spectrum of the total emission (pulsed+unpulsed) and 
phase-resolved spectra of the pulsed emission. We use 
five equally spaced phase intervals. NuSTAR data are 
fitted above 10 keV, where the coronal outflow has to 
account for most of the observed emission. 

Figure [TUI shows the map of the p- value of the fit in the 
a m ag-/3obs plane. The acceptable model is clearly iden¬ 
tified in this map0 It has a mag = 0.25 ± 0.15, /3 0 bs = 

4 There are in fact two solutions because interchanging the values 
of a m ag and /3 0 b s does not change the model spectrum, as long as 
the j-bundle is assumed to be axisymmetric. 


1.0 ± 0.2, 9j = 0.24 ± 0.02, and A 9j/9j > 0.26, consis¬ 
tent with a polar cap. The corresponding magnetic flux 
in the j-bundle is (2.5 ± 0.4) x 10 26 G cm 2 . The power 
dissipated in the j-bundle is Lj = (6 ± 1) x 10 36 erg s _1 . 
Most of the released energy is radiated in the MeV band 
(peaking at ^6 MeV) and is not seen to NuSTAR. Us¬ 
ing the obtained best-fit model for the hard X-ray com¬ 
ponent, we have investigated the remaining soft X-ray 
component. The procedure is similar to that in A13 and 
lHascoet. Beloborodov fc den Hartog I (2014); we freeze 
the best-fit parameters of the outflow model and fit 
the spectrum in the 0.5-79 keV band using the NuS¬ 
TAR., Swift, Chandra, and XMM-Newton data. As in 
A13, we find that the 0.5-79 keV spectrum is well fit¬ 
ted by the sum of two blackbodies (which dominate be¬ 
low 10 keV) and the coronal outflow emission (which 
dominates above 10 keV). The cold and hot blackbod¬ 
ies have luminosities L c = 2.2 ± 0.1 x 10 35 erg s _1 , 
Lh = 9.8 ± 1.3 x 10 34 erg s _1 and temperatures kT c = 
0.45 ±0.01 keV, kT h = 0.75 ±0.02 keV. Note that these 
values are different from those we obtained with the phe¬ 
nomenological models in Table 0] 
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Figure 10. Map of p-values for the fit of the hard X-ray compo¬ 
nent with the coronal outflow model; the p-values are shown in the 
plane of (u mag , /3 0 b s ) anil maximized over the other parameters. 
The p-value scale is shown on the left. The hatched green re¬ 
gion has p-values smaller than 0.001; the white region has p-values 
greater than 0.1. Interchanging the values of « m ag and 0 o \, a does 
not change the model spectrum, as long as the j-bundle is assumed 
to be axisymmetric. Therefore, the map of p-values is symmetric 
about the line of /3 0 b s — a mag- 

4. SWIFT MONITORING OBSERVATIONS 

We report below on Swift monitoring observations for 
spectral and temporal behavior of the source on long 
timescales. The observations were taken with the X-Ray 
Telescope (XRT) using Windowed-Timing (WT) mode 
for all observations, which have been conducted once ev¬ 
ery two to three weeks since 2011 July, except when the 
source was in Sun-constraint from mid-November to mid- 
February each year. In total, 68 observations (not listed 
in Table 1) having ~266 ks of summed exposure were an¬ 
alyzed. The Swift data were processed with xrtpipeline 
using the HEASARC remote CALDB0 


4.1. Data Analysis and Results for the Swift Monitoring 
Observations 

To investigate the spectrum of IE 1841—045 in the 
monitoring observations, we extract spectra for each ob¬ 
servation using a 10-pixel (24") long strip centered on 
the source. An annulus of inner radius 75-pixel, and 
outer radius 125-pixel centered on the source was used to 
extract background spectra. The spectra were grouped 
to have a minimum of 20 counts per bin. The spec¬ 
tra were fitted with a photoelectrically absorbed black- 
body with an added power-law component, using the 
tbabs (bbody+pow) model in XSPEC 12.8.1, with Ah held 
fixed at 2.05 x 10 22 cm 2 , the value we obtained in Sec¬ 
tion [3223 No significant change in source 1-10 keV flux 
was observed over the monitoring period of ^3 years 
(x 2 /dof=49/59), including the AmSTA ill-observed burst¬ 
ing period. The result is shown in Figure [TT] (a). 

We also searched all the Swift observations for bursts 
by binning the source region light curves above 1 keV 
into 0.01 s, 0.1 s, and 1.0 s bins. The counts in each 
bin were compared to the mean count rate of its GTI, 
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Figure 11. Results of the Swift monitoring campaign for flux 
and timing, (a) 1-10 keV flux, (b) timing residuals after fitting 
out v and z>, and (c) after fitting out twelve frequency derivatives. 
Vertical dashed lines show the periods when NuSTAR observations 
were taken. 

assuming the Poisson distribution. We found no signifi¬ 
cant bursts in the Swift XRT data. Note that the Swift 
observations did not cover the NuSTAR- detected bursts 
times presented in Table [2] 

In order to derive the timing solution and to search 
for glitching activity, we barycentered the source events 
using the location of IE 1841-045, R.A.=18 ,l 41 m 19.34 s , 
Decl.=—4°56 , ll , /2. We then extracted Time-Of-Arrivals 
(TOA) using a Maximum Likelihood (ML) method (see 
iLivingstone et al.l 120091 iScholz et al.l I2012T ) . The ML 
method compares a continuous model of the pulse pro¬ 
file to the profile obtained by folding a single observa¬ 
tion. These TOAs were fitted to a pulse arrival model 
(e.g., the q uadratic fu ncti on i n Section 13.2. ID using the 
TEMP02 (jHobbs. Edwards fe Manchester! 120061) pulsar 
timing software package. We find that a timing model 
consisting of v and v does not fit the data well (Fig¬ 
ure [TT] (b) and “Fit 1” in Table [5]) as was already ob¬ 
served in the previous RX TE monitoring observations 
(e.g., see lDib fc Kaspil[2014f) . We need to include twelve 
frequency derivatives in order to achieve an acceptable fit 
(i.e., x 2 /dof~l) with Gaussian residuals (Figure [Tl] (c) 
and “Fit 2” in Table [5]). Note that the timing solutions 
presented in Table [5] are valid only over the time interval 
of the monitoring campaign. 

Motivated by the apparent ‘kink’ in the residuals of 
the simple spin-down model around MJD 56100, we at¬ 
tempted to fit a glitch at the epoch of the kink. However, 
the data are better fit using a model with four frequency 
derivatives (rms residual of 0.97 s) versus a glitch model 
(rms residual of 0.99 s). Therefore we do not need to 
invoke a sudden event to explain the measured TOAs. 
We present our best timing solutions in Table [5] 


1 http: //heasarc.nasa.gov/docs/heasarc / caldb / caldb_remote_acce 
ss.html 5. DISCUSSION 
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Table 5 

Timing Parameters for IE 1841—045. 


RA (J2000) 

18 h 41 m 19.34 s 

DEC (J2000) 

—4°56'11.2" 

MJD Range 

55795-56799 

Epoch (MJD) 

56300 

Fit 1 

* (s’ 1 ) 

0.084 806 860 6(9) 

+ 

CO 

1 

to 

—2.9121 (8) x 10“ 13 

RMS (s) 

4.82 

X 2 /dof 

2889.05/53 

Fit 2 

* (s' 1 ) 

0.084 806 897(4)) 

CO 

1 

to 

—2.985(12) x 10“ 13 

# (s- 3 ) 

6.8(23) xlO- 22 

# (s’ 4 ) 

4.6(8)xl0 -28 

& ( s “ 5 ) 

—8.7(16)xl0 -35 

# (s- 6 ) 

—2.5(5)xl0 -41 

# (s- 7 ) 

6.0(12) xlO -48 

# (s- 8 ) 

1.2(3) xl0“ 54 

^ (s- 9 ) 

-3.4(7) xlO" 61 

# (s- 10 ) 

-4.5(10) xlO -68 

O’- 11 ) 

1.4(3)xl0 -74 

(s- 12 ) 

9.5(22) xl0~ 82 

& (s- 13 ) 

—3.3(8)xl0 -88 

RMS (s) 

0.43 

X 2 /dof 

49.34/42 


Notes. All errors are TEMP02 reported lcr errors. 

The new 350-ks observation of IE 1841045 by NuS¬ 
TAR allowed a significantly better study of its persistent 
emission and the serendipitous detection of X-ray bursts. 
Below we discuss the results and compare them with ob¬ 
servations of other magnetars. 

5.1. The X-ray Bursts and the Tails 

Magnetars often show bursting behavior in the X-ray 
band, which may be caused by i nstabi li ties inside the 
neutron star or i ts magnetosphe re ([Thompson fe Duncaiil 
119951 ILvutikovI 120031 Woods fe Thompsonl I2006lh The 
time pro files and spec t ra of bursts show significant di¬ 
versity. IWoods et al.l (1200511 suggested that there are 
two types of magnetar bursts, one having significant 
tail emission and the other having orders of magnitude 
smaller tail emission. The authors attributed the for¬ 
mer to cr ustal activ ity a nd the latter to magnetospheric 
activity. iLenters et al.1 (12003T ) found a strong correla¬ 
tion between bursts and tail energies in the magnetar 
SGR 1900+14. 

The X-ray bursts from IE 1841—045 in 2013 Septem¬ 
ber have very short rise and fall times, with Xgo of 
0.01-0.6 s, and hard spectra (T=l-2 or kT ~3-5 keV). 
The blackbody temperature we measured for the burst 
5 (4.8 + 0.5 keV; see Table [3]) is consistent with that of 
the colder bla ckbody measured with the GBM (fcT j = 
5.3 ± 0.2 keV; [Collazzi . Xiong, fe Kouveliotoull2014D . 

iKumar fc Sah-Harbl (|2010D reported detection of emis¬ 


sion lines at 27 keV, 40 keV and 60 keV in the 2010 May 
burst spectrum with Swift BA T, alt hough these are ar¬ 
gued against later bv lLin et al.l (|2011l) . Interestingly, all 
the lines are in the NuSTAR band, and could be detected 
by NuSTAR if they appeared again. However, we do not 
see evidence of line emission. Therefore, we estimate 
90% upper limits on any 27keV Gaussian line flux to be 
0.24 photons cm~ 2 s -1 and 0.11 photons cm~ 2 s _1 in 
the brightest burst spectrum for the blackbody and the 
power-law continuum models, respectively. Note that 
we are not able to comp are our results with those of 
IKumar fc Sah-Harbl (2 0101 ) since they did not present the 
line flux. 

An extended tail is reliably detected only in ener¬ 
getic burst 5, and we find a hint of a tail in burst 6. 
Note that these two bursts were also observed by the 
GBM, and had si g nificant flux a bove t he NuSTAR band 
(ICollazzi. Xiong, fe Kouveliotoull2014D : they are the two 
most energetic bursts in our sample. Thus, the tail 
brightness and the burst energy we measure for the 
IE 1841—045 bursts seem to agree qualitatively with the 
correlation reported for SGR 1900+14 a nd SGR 1806—20 
(jLenters et alJl2003L iGogiis et al.l 1201 111 . A similar trend 
was also seen in the recent bursts from IE 1048.1—5937 
(lAn et al.ll2014f) . We note however that the large energy 
seen by GBM in bursts 5 and 6 indicates a tail-to-burst 
ratio (-Stall/Eburst ~ 2 x 10 -2 ) that is much lower than in 
IE 1048.1—5937 (-Etaii/Eburst ~ 5-60). This confirms the 
known diversity of tails o f magnetar bursts (|Kasni et al.1 
12004 IWoods et, aD 12005th Whether the tail-to-burst ra¬ 
tio follows a bimodal or a random distribution is not yet 
clear, and further investigation is required. 

In contrast to the burst tails observed by NuSTAR for 
IE 1048.1—5937, the burst tail in IE 1841—045 shows 
no clear correlation between spectral hardness and flux. 
This correlation was also absent in some of the long¬ 
term (mo nths to years) f lux relaxation of other magne¬ 
tars (e.g., lAn et al.l[20m 

The flux evolution in the tail of burst 5 followed a 
power-law decay with a decay index of 0.45 ± 0.10. 
The flux decay is similar to the tails o f bursts from 
SGR 1900+14 (decay index of 0.43-0.7; ILenters et all 
120031) . A significantly faster decay was observe d for burst 
tails in IE 1048.1—5937 (decay index of 0.8-1: lAn et al.1 
12014 . It is possible that the tail contains many unre¬ 
solved weaker bursts which affect the observed flux decay. 
It is still unclear what controls the resulting decay rate 
and why it is significantly different in IE 1048.1—5937 
and IE 1841-045. 

Furthermore, we find that the tail spectra in 
IE 1841—045 are similar to (or slightly softer than) 
the persistent emission. In contrast, the tail spectra in 
IE 1048.1—5937 were significantly harder than its per¬ 
sistent emission. This further contributes to the diver¬ 
sity of magnetar bursts. For instance, IE 2259+586 ex¬ 
hibited bursts with and without ks-long tails, and the 
tai l emission was o bserved to soften with decreasing flux 
(jKaspi et al.1 12004 ). The observed diversity of magnetar 
bursts is not well explained by current theoretical mod¬ 
els. 

5.2. Pulse Profile and Pulsed Fraction 

It is known that the pulse profiles of magnetars 
can look significantly different in different energy 
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bands (jden Hartog et al.l 12008ft . This is also true for 
IE 1841—045 (Figure [5]). An interesting feature of the 
pulse profile is the double-peaked structure in the nar¬ 
row band of ~24-35 keV. A13 found this in the previ¬ 
ous NuSTAR observation, and suggested that it may be 
caused by a spectral feature. The new long observation 
confirms the double peaked structure. It shows however 
that the change in the pulse profile is not as sharp as was 
suggested previously, and its shape may be more compli¬ 
cated, as a hint of another peak appearing between the 
two peaks is seen in several energy bands (e.g., in the 
33-35 keV profile in Figure 0 . 

The pulsed fraction (and the rms pulse amplitude) in¬ 
creases with photon energy below ~8 keV (Figure 0 . 
We found no significant increase in the pulsed fraction 
above 8 k eV, i n co ntrast with previous measurements 
(jKuiner et al.l 12006 1. This conclusion is insensitive to 
the choice of energy bins in Figure 0 which can affect 
the measured values in individual bins, but not the gen¬ 
eral trend. Our spectral analysis also suggests that the 
pulsed and steady components of the flux have similar 
power-l aw p hoton indices at high energies >15 keV (see 
Section 13.31 and Tabic 0 . In the coronal outflow model, 
this implies that one of the angles, a mag and /3 0 bs, is 
small, or the emission region ( 0j ) is broad. 

5.3. Spectra 

We have measured the phase-averaged, phase-resolved 
and pulsed spectra of the source. Our best-fit parameters 
are slightly different from those reported by A13. The 
use of the new CALDB for the NuSTAR data may have 
some effect on the obtained spectral shape, however we 
showed in Section f3. 2. 31 that this effect is small. 

We do not see any significant change in the source 
flux among the observations taken over one year de¬ 
spite the fact that bursts were detected in some obser¬ 
vations but not in the others. We further compared the 
NuSTAR-me&sured spectrum with that reported previ¬ 
ously (r = 1.32±0.11 and Tm-inokfi V = 3.0 x 10 35 erg s" 1 
for an assumed distance of 6.7 kpc; iKuiper et al.ll2006ft . 
and find that our measurement (T = 1.37 ±0.01 and 
Tio-iookeV = 3.0 x 10 35 erg s -1 ) is fully consistent with 
the previous values, suggesting that the source hard-band 
spectrum has been stable over 10 years. The same trend 
in the soft band (1-10 keV) is seen in the Swift mon¬ 
itoring data (see Figure [FT]). Stability in the soft- band 
pulsed flux has been reported bv lDib fc Kasnil (|2014ft and 
Zhu fc Kasnil (12010ft for a longer period. We note that 
Kumar fc Safi-Harbl (i2010ft reported a possible increase 
in the soft-band flux for IE 1841—045 (<10 ke V) associ¬ 
ated w ith the source burst activity. However, iLin et ahl 
(1 20 Il f) showed that there is no significant change in the 
soft-band flux over a 1 400-day i nterval , incl uding the 
burst period reported bv lKumar fc Safi-Harbl 1 2010 '). 

The new measurements of the phase-resolved spectrum 
agree with the previous NuSTAR observation. The spec¬ 
trum is harder near the pulse peaks, with smaller Th and 
greater Ebreak- We also note a possible hardening of the 
phase-averaged spectrum above >15 keV. However, we 
cannot reliably measure the spectrum curvature in the 
hard X-ray band due to the statistical noise. 

We have searched for the spectral feature that was 
suggested as a possible explanation for the change in 
the pulse profile near 30 keV (A13). We find that the 


phase-resolved spectra are all well fit by a blackbody 
plus broken power law or a blackbody plus double power 
law; no emission or absorption line is required or hinted 
at by the fit residuals. We also searched for a spectral 
feature that shifts with rotation al phase as was seen in 
SGR 0418±5729 dTiengo et al.ll2013 i. but did not find 
such a feature (Figure 0. 

5.4. Outflow model 

The more accurate phase-resolved spectrum obtained 
from the 350-ks NuSTAR observation provides a new op¬ 
portunity to apply the coronal outflow model. We found 
that the model still fits the observed hard X-ray spec¬ 
trum, and it does so in a small region of the parameter 
space (Figure fTT)l). which allowed us to estimate the size 
of the active j-bundle, the dissipated power in the mag¬ 
netosphere, and the angles between the magnetic axis, 
the rotation axis, and the line of sight. The increase in 
exposure by a factor ^6 excluded at more than 3cr level 
the second solution for the outflow model that was found 
in A13. 

Interestingly, adding the new free parameter A 0j does 
not introduce new acceptable solutions (with a p-value 
above 10 -3 ) and does not significantly affect the best 
solution — the best fit shows that the footprint of the 
j-bundle on the star is a broad ring, hardly distinguish¬ 
able from a polar cap. This contrasts with IE 2259±586, 
where a thin ring with A 9j/6j < 0.2 is clearly statisti¬ 
cally favored over a polar cap. This diversity may point 
to different distributions of the crustal magnetic stresses 
which are responsible for the magnetospheric twisting. 

The soft X-ray component (below ~10 keV) is well fit¬ 
ted by the sum of two blackbodies. The best-fit model is 
similar to that found in A13. The cold blackbody covers 
a large fraction of the star area, A c ~ 0.42AnS: where 
Ans is the area of the neutron star with an assumed ra¬ 
dius Ans = 10 km. The emission area of the hot black¬ 
body, Ah ~ 0.024>Ins, is comparable with the area of 
the outflow footprint Aj = 7rsin 2 9j ~ 0.014 (6j/ 0.24) 2 . 
This is consistent with the coronal outflow model, where 
the footprint of the j-bundle is expected to form a hot 
spot, as some particles accelerated in the j-bundle flow 
back to the neutron star and bombard its surface. 

Figure [TO] suggests that the coronal outflow correctly 
describes the hard X-ray source as a decelerating 
outflow ejected from a discharge zone near the star. 
The outflow parameters inferred from the fit of the 
phase-resolved data have rather small statistical uncer¬ 
tainties, and the results reveal a puzzli ng feature of 
IE 1841—045. Using Equation (48) in IBeloborodovI 
(12009ft . one can infer the discharge voltage in the active 
j-bundle: $ ~ lO 11 ^ -1 V where 'P is the twist implanted 
i n the magnetosp h ere, w hich does not exceed V'max ~ 3 
dParfrev. Beloborodov fc Huill2013i j. Due to our refined 
constraint on 0j and the strong dependence of $ on 8j 
(4> oc 9j ), the inferred voltage is one order of magni¬ 
tude higher than that given in A13. The high voltage 
is surprising in two ways: (1) it exceeds the expected 
thres hold for e ± discharge ((Beloborodov fc Thompson! 
12007ft by at least a factor of 10; (2) it implies a short 
timescale for ohmic dissipation of the m agnetospheric 
twist fdiss ~ 0.1 ip 2 yr (Equation (50) in IBeloborodovI 
12009ft . Thus, without continued energy supply, the j- 
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bundle is expected to untwist on a year timescale or 
faster, which is not observed — the persistent X-ray 
emission from IE 1841—045 has been stable for at least 
one decade. 

This puzzle is related to a more general question: why 
are some magnetars transient and others persistent? The 
magnetospheric activity should be fed by magnetic en¬ 
ergy pumped from the star by its surface motions. The 
surface motions are caused by the crust yielding to ac¬ 
cumulating magnetic stresses inside the magnetar. En¬ 
ergy supply to the twisted magnetosphere may or may 
not be intermitte nt, depend ing on the mechanism of the 
crustal motions. iBeloborodov fc Levin I (120141) have re¬ 
cently shown that the crust can yield through a thermo¬ 
plastic instability which launches a slow wave resembling 
the deflagration front in combustion physics. The ther¬ 
moplastic wave rotates the crust and twists the magneto¬ 
sphere. This mechanism naturally triggers the outbursts 
observed in transient magnetars, yet it is still unclear 
how the quasi-steady state observed in IE 1841—045 is 
formed. It could be formed by frequently repeated en¬ 
ergy supply to the j-bundle combined with some form of 
feedback on its untwisting rate. 

5.5. Swift monitoring observations 

We did not detect any changes in the 1-10 keV source 
flux during the three years of Swift-X RT monitoring since 
2011. In particular, the flux did not change significantly 
during the bursting period, which is consistent with the 
stability of persistent emission in the NuSTAR data. Our 
measurements do not agree with the observation that 
the sou rce persistent emissio n properties chan ged due to 
bursts dKurnar fc Safi-Harbl 12010 ). Note that Lin et all 
( 2011 1 also observed no change in the persistent proper¬ 
ties of the source after bursts. The flux stability also im¬ 
plies that the source seems to be in a perpetually burst¬ 
ing state, and perhaps this is a difference between the 
classical bright AXPs monitored with RXTE and ‘tran¬ 
sient’ sources, i.e., perhaps the bursting indicates that 
there is more heat dissipating and hence more magnetic 
activity in this magnetar, consistent with the higher Lx 
in “quiescence”. 

The pulse t iming behavior of IE 1841—045 is known to 
be noisy, and lDib fe Kaspil (|2014H had to use five deriva¬ 
tive terms to fit the timing data. Our timing solution re¬ 
quires twelve freque ncy derivatives. Note that although 
iDib fc Kaspil ( 2014 1 used a much longer data span (13 
yr), they break the data into smaller pieces which span 
1-3 yr, similar to ours. The Swift monitoring data re¬ 
quire more frequency derivatives than the RXTE data 
do due to the large kink at MJD 56100 which may in¬ 
dicate a glitch. We investigated the possibility that the 
kink is caused by a glitching activity but did not find 
any evidence of a glitch during the monitoring period of 
3 years including the bursting period. 

6. CONCLUSIONS 

During NuSTAR observations of the magnetar 
IE 1841—045, we detected six X-ray bursts from the 
source. The bursts are short, Tgo < 1 s, and bright. 
A tail was observed after one burst which was most en¬ 
ergetic as measured with Fermi GBM. The tail emission 
was similar to or softer than the persistent emission, and 


the flux decay in the tail followed a power law with a de¬ 
cay index of ^0.5, with no clear spectral softening with 
time. The properties of the tail emission are different 
from those seen after recent NuSTAR-observed X-ray 
bursts from the magnetar IE 1048.1—5937 whose tails 
decayed fast with flux decay indices of 0.8-0.9 and had 
harder spectra than the persistent emission. The new 
observations also yield detailed pulse profiles in different 
energy bands, and show that the pulsed fraction does not 
increase rapidly with energy, in contrast to previous ob¬ 
servational reports. We show that the source hard-band 
flux has been stable over ^10 years. Using a Swift mon¬ 
itoring campaign, we found that the 1 10 keV flux from 
the source has been stable within <20% and the source 
timing behavior has been very noisy during a 3-year pe¬ 
riod, in spite of the bursting behavior we have observed. 

The X-ray spectra of IE 1841—045 are well fitted by 
the coronal outflow model of IBeloborodov! (2013 '). The 
new fit has provided improved constraints on the angle 
between the magnetic and rotation axes of the magnetar 
and the size of its active j-bundle. Remarkably, the best 
fit implies fast dissipation of the magnetic twist, in ap¬ 
parent contradiction with the observed stability of X-ray 
emission, as no long-term evolution has been detected in 
the persiste nt emission of IE 1 841—045 for more than 
one decade (iKuiper et al.ll2006l A13). This behavior is 
distinct from the untwisting magnetospheres in transient 
magnetars and presents a puzzle that is yet to be re¬ 
solved. 
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APPENDIX 

COMPARISON BETWEEN ESTIMATORS OF PULSED FRACTION 


Some confusions on pulsed fraction of a pulsar have been arisen mainly because different literatures use a different 
estimator. In this appendix, we review several commonly used estimators of pulsed fraction, and describe cautions to 
be taken when using them. 

We consider four different estimators of pulsed fraction commonly used in literatures: pulsed fraction measured by 
(1) count area PF area (Equation [2]), (2) min-max counts, PF m i nmax = Pmax ~ Pmin , (3) fitting the light curve and 

PmaxT Pmin 


(4) calculating rms variation 




12i(Pi-<P>) 2 ' 

<p > a i 


using truncated Fourier series and subtracting Fourier power produced 


by noise in the pulse profile, a a k +&l k , from the Fourier amplitudes (PF r ms , Equation^ see also I Archibald et ah]l2014lh 
As we show below, estimators (1) and (2) are biased upwards since one has to select the minimum (and maximum) 
count bin in the pulse profile. For example, if the DC component of a pulse profile is relatively broad, one of the DC 
bins will very likely fall below the “true” minimum because of statistical fluctuation. Since one will use that phase 
bin for p m i n , the resulting pulsed fraction will be higher than the true value. Although this bias can be avoided by 
holding the minimum phase ((/> m in) fixed and performing simulations, the spread in the measurements will be larger 
in this case than in the case of picking up the minimum bin as we show below. 

We conducted simulations to compare the estimators. For the simulations, we used a sine plus constant function 
(Asin^> + B) for the theoretical pulse profile. The pulsed fraction for this profile can be analytically calculated for all 


the estimators, and is A/B for (1)—(3), and A/(Bv / 2) for (4). We carried out 10,000 realizations of the pulse profile 
for A = 25 and B = 50 per phase bin having a total of ^800 events in the pulse profile (Figure [T3] (a)), and measured 
the pulsed fraction using the above estimators. 

After each realization, we measured the pulsed fraction as defined above. Note that for the estimator (3), we used 
a sine plus a constant function for the fit. We show the results in Figure [T3] and Table [G] We find that estimators 
(1) and (2) give a biased value with a chance of 98% and 99.8%, respectively (see Figures [TH (b) and (c)). When 
holding the minimum phase fixed at the theoretical minimum phase, no significant bias is seen for both the estimators, 
but the spread of the measurements increased (Tabic [6]). The other estimators, (3) and (4), provide unbiased results 







17 



Figure 12. Simulation results: (a) pulse profiles where the solid line is a theoretical pulse profile, and the dotted curve shows a realization 
of the theoretical pulse profile, and pulsed fractions measured using (b) area fraction (PF are a), (c) min-max bin (PF m i nmax ), (d) a function 
fit (PFfit), (e) rms amplitude (\/2PFrms), and (f) for different binning. Vertical red dotted line in panels (b)-(f) shows the theoretical 
value for the pulsed fraction. 


Table 6 

Summary of measurements made with 10,000 simulations 


Estimator 

Mean 

Spread (1 <t) 

A a 

^theory 

Notes 

PFarea 

0.61 

0.06 

0.11 


P Farea 

0.5 

0.1 

0.0 

Fixed 0 m i„ 

PF • 

1 1 minmax 

0.63 

0.05 

0.13 


PF 

1 1 minmax 

0.50 

0.09 

0.0 

Fixed 0 min 

PF { it 

0.501 

0.031 

0.001 


P Frms 

0.512 

0.032 

0.012 


P Frms 

0.513 

0.032 

0.013 

No correction term 


a Difference between the theoretical value and the mean of the measurements. 


(Figures 0 (d) and (e)). We also measured PF lms in equation [3] without the small bias-correction term, a^ k + a kk 
which is to remove the underlying Fourier power produced by noise in the pulse profile. We find that this term changes 
the results only by ~0.1% (see Table[6]). 

We carried out simulations with different parameters (A and B ) for the theoretical pulse profile and find that the 
results change depending on the number of events and the pulse fraction. For example, the estimators (1) and (2) 
tend to provide better results as the total count or the pulsed fraction increases. However, the upward bias does 
not disappear even for a fairly large number of total counts (e.g., 10,000) or pulsed fraction (e.g., 70%). The other 
estimators did not significantly bias the results in any set of parameters we studied. 

We investigate effects of binning as the results can change depending on binning. We binned the light curve into 8, 
16, 32, 64, 100, and 128 bins. Note that we changed the total number of events for different binning to keep the average 
number of counts per bin same and to have similar statistical error for the minimum (and maximum) phase bin. The 
results are shown in Figure fl2l (f). We find that estimators (1) and (2) bias the results larger as the number of bins 
increases while (3) and (4) produce robust results regardless of binning. This is expected as there are more bins among 
which we can choose the minimum when we increase the number of bins, and hence it is more likely for estimators 
(1) and (2) to have p m in smaller than the true minimum. Estimators (3) and (4) do not rely on one minimum phase 
bin but on the statistical average of the DC component, and thus they are insensitive to the number of bins used. 
Furthermore, they provide more accurate results as the number of bins increases simply because we simulated more 
events in the cases of finer binning. 

We find that estimator (3) provides the best results; the measurements are closest to the theoretical value, and the 
spread is smallest. Note that we had a priori knowledge on the pulse shape with which we fit the pulse profile for 
that estimator. Furthermore, the simple pulse shape allows us to fit the profile only with two parameters. If the pulse 
shape were more complex, one may have to use more harmonics to fit the pulse profile and the results may be similar 
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to those obtained using PF lms (e.g., see Figure [Gl for error bars in two-harmonic fits). Hence, it may be easier to use 
PFims estimator for more complex pulse profiles, although the value measured with this estimator is different from 
that of the others in general, and the conversion factor (e.g., \/2 for the sine function in the above simulation) can 
change depending on the pulse shape, making direct comparison with the others difficult. 

In summary, we find that the pulsed fraction estimators PF area and PF r ninmax, often used in literatures, give a 
biased result in general and are sensitive to the number of bins used. Therefore, they should be used with great care. 
The other estimators, PFa t and PF lms , provide an accurate measurement regardless of binning, and hence they are 
preferred. We note, however, that results of PF rms cannot be directly compared with those of the others since there 
is a scale factor which differs for different pulse shape; the factor can be obtained using simulations. 

One final remark we would like to make is that extra care needs to be taken when comparing pulsed fractions 
measured with different instruments. In particular, if the pulse fraction changes strongly over the energy band it was 
measured, the result should be regarded as an energy-weighted pulsed fraction. The energy-weighted pulsed fraction 
will be measured differently by different instruments since they have different energy responses. In this case, one has 
to use either a smaller energy range over which the pulsed fraction does not change much or a response-unfolded 
estimator such as flux density ratio (A13). 


